Insulating Behavior of an Amorphous Graphene Membrane 
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We investigate the charge transport properties of planar amorphous graphene that is fully topologically dis- 
ordered, in the form of sp 2 three-fold coordinated networks consisting of hexagonal rings, but also including 
many pentagons and heptagons distributed in a random fashion. Using the Kubo transport methodology and 
the Lanczos method, the density of states, mean free paths and semiclassical conductivities of such amorphous 
graphene membranes are computed. Despite a large increase in the density of states close to the charge neutral- 
ity point, all electronic properties are dramatically degraded, evidencing an Anderson insulating state caused by 
topological disorder alone. These results are supported by Landauer-Biittiker conductance calculations, which 
show a localization length as short as 5 nanometers. 

PACS numbers: 72.80. Vp, 73.63.-b, 73.22.Pr, 72.15.Lh, 61.48.Gh 



The physics of disordered graphene is at the heart 
of many fascinating properties such as Klein tunnel- 
ing, weak antilocalization or anomalous quantum Hall ef- 
fect (see reviews 1 ^). The precise understanding of in- 
dividual defects on electronic and transport properties of 
graphene is currently of great interest^. For instance, 
graphene samples obtained by large-scale production meth- 
ods display a huge quantity of structural imperfections 
and defects that jeopardize the robustness of the other- 
wise exceptionally high charge mobilities of their pristine 
counterparts 4 .. Indeed, the lattice mismatch-induced strain 
between graphene and the underlying substrate generates 
polycrystalline graphene with grain boundaries that strongly 
impact transport properties^. However, despite the large 
amount of disorder, such graphene flakes (when deposited 
onto oxide substrates) usually maintain a finite conductiv- 
ity down to very low temperatures owing to electron-hole 
puddles-induced percolation effects that preclude localiza- 
tion phenomena close to the Dirac point^. The predicted An- 
derson localization in two-dimensional disordered graphene 
has been hard to measure in non intentionally damaged 
graphene, in contrast to chemically modified graphene 7 ^. 
Nevertheless, in a recent experiment it was possible to screen 
out electron-holes puddles using sandwiched graphene in be- 
tween two boron-nitride layers, together with an additional 
graphene control layer 9 . As a result of puddles screen- 
ing, a large increase of the resistivity was obtained at the 
Dirac point, evidencing an onset of the Anderson localiza- 
tion regime. 

Beyond individual defects and polycrystallinity, a higher 
level of disorder can be induced on graphene to the point of 
obtaining two-dimensional fully amorphous networks com- 
posed of sp 2 hybridized carbon atoms. Such networks con- 
tain rings other than hexagons in a disordered arrangement. 
The average ring size is six according to Euler's theorem, 
allowing these systems to exist as flat 2D structures. Ex- 
perimentally, such amorphous two-dimensional lattices have 



been obtained in electron-beam irradiation experiment s 10 - 11 , 
and directly visualized by high resolution electron trans- 
mission microscopy. Previously, indirect evidence for the 
formation of an amorphous network was obtained using 
Raman spectroscopy in samples subject to electron-beam 
irradiation' 2 , ozone exposure^ and ion irradiation^ 4 .. In 
all these cases, an evolution from polycrystalline to amor- 
phous structures was observed upon increase of the dam- 
age treatment. Inii, further evidence of the formation of 
an amorphous network was obtained through transport mea- 
surements. These indicate the transition from a weak lo- 
calization regime in the polycrystalline samples to variable 
range hopping transport in the strongly localized regime for 
amorphous samples, as evidenced by the temperature depen- 
dence of the conductivity. Localization lengths were esti- 
mated to be of the range 0.1 to 10 nm in the amorphous 
samples, depending on the degree of amorphization. From 
the theoretical side, models of the amorphous network have 
been proposed using stochastic quenching methods 1 ^, and 
molecular dynamics 1 ^—. Electronic structure calculations 
show that the amorphization yields a large increase of the 
density of states in the close vicinity of the charge neutrality 
pointer—. Despite the expected reduction of the conduc- 
tion properties due to strong localization effects, Holmstrom 
et alj& suggest that disorder could enhance metallicity in 
amorphized samples, in contrast with the experimental evi- 
dence. 

In this Letter, we explore the transport properties of two- 
dimensional sp 2 lattices with a massive amount of topolog- 
ical disorder, encoded in a geometrical mixture of hexagons 
with pentagon and heptagon rings. The calculations are done 
using two complementary approaches: a Kubo formulation 
in which the conductivity of bulk 2D amorphous graphene 
lattices is determined, and a Landauer-Biittiker formulation 
where the conductance of stripes of amorphous graphene 
contacted to semi-infinite pristine graphene electrodes is cal- 
culated. Both approaches lead to similar findings. Depend- 
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ing on the ratio between odd versus even-membered rings, a 
transition form a graphene-like electronic structure to a to- 
tally amorphous and smooth electronic distribution of states 
is obtained. The stronger the departure from the pristine 
graphene, the more insulating is the corresponding lattice, 
which transforms into a strong Anderson insulator with elas- 
tic mean free paths below one nanometer and localization 
lengths below 10 nm close to the charge neutrality point. 
Those structures are therefore inefficient to carry any sizable 
current, and unsuitable for practical electronic applications 
such as touch screens displays or conducting electrodes, but 
interesting for scrutinizing localization phenomena in low 
dimensional materials. 



TABLE I: Characteristic parameters for the two samples of amor- 
phous graphene, SI and S2 




FIG. 1: (color online), (a), (b) Details of amorphous graphene sam- 
ples S 1 and S2, respectively, used to compute the conductivity with 
the Kubo approach, (b) Scheme of setup for the Landauer-Biitikker 
calculations. The grey area represents the amorphous graphene 
stripe. Periodic boundary conditions are used in the y direction. 
The graphene electrodes are semi-infinite in the x direction, (d) 
Total density of states of the two amorphous samples. The pristine 
crystalline graphene case (dashed lines) is also shown for compar- 
ison. 



Models of amorphous graphene.- Amorphous models 
of graphene are prepared using the Wooten- Winer- Weaire 
(WWW) methodic, introducing Stone-Wales defects 2 ! 
into the perfect honeycomb lattice. These networks can 
be considered as the amorphous versions of the Haeckelite 
structures proposed by Terrones et al.—. To generate the 
structures, periodic boundary conditions are imposed and the 
entire network is relaxed with the Keating-like potentia l 15 ' 23 . 
A further relaxation of these models using forces from Den- 
sity Functional Theory leads to slight changes in the bonding 
distances and angles, but to virtually identical radial distribu- 
tions functions and, most importantly, electronic properties 





SI 


S2 


number of atoms 


10032 


101640 


% of n-membered rings (n = 5/6/7) 


24/52/24 


44/12/44 


< n 2 > - < n > 2 


0.47 


0.88 


RMS deviation of bond angles 


11.02° 


18.09° 


RMS deviation of bond lengths 


0.044 A 


0.060 A 


Fermi energy (in units of 70) 


0.03 


0.05 



like the density of states, as shown in Ref. [l7l Fig. Qla)-(b) 
show pieces of the two studied sample S 1 and S2, which re- 
spectively contain 10032 and 101640 atoms (all of them with 
three-fold coordination as the honeycomb lattice, but topo- 
logically distinct). The amorphous character of the samples 
is demonstrated byanalyzing the radial distribution function, 
as shown in Ref. 15} Table|I]shows the parameters that char- 
acterize the two samples. For sample S 1 , 24% of the elemen- 
tary rings are pentagons, 52% hexagons and 24% heptagons, 
while sample S2 has a larger share of odd-membered rings. 
In both samples, the number of heptagons is the same as that 
of pentagons, as required by Euler's theorem, and these sys- 
tems can exist without an overall curvature as flat 2D struc- 
tures with some distortions of bond lengths and angles, al- 
though they may pucker under some circumstances. We will 
only be concerned with the planar structures here. Sample 
S2 is an extreme case, having very few hexagons and being 
furthest from the pristine honeycomb lattice. This is useful 
to accentuate the differences between crystalline and amor- 
phous samples, and to gain perspective. Nevertheless, it is 
likely that sample SI is nearer to physical reality as it is less 
strained. 

For the calculation of the Landauer-Biittiker conductance, 
we set up models in which an amorphous stripe is contacted 
by two pristine graphene electrodes separated by a distance 
L, as shown in Fig. QI C )- We use models with increasing 
values of L, to explore the dependence of the conductance 
on the length of the amorphous contact in the transport di- 
rection. The models are periodic in the direction perpendic- 
ular to the stripe, with a periodicity of W=\ 1.4 nm, and they 
have the same ring statistics as the bulk sample S 1 described 
above. 

Electronic Properties .- The electronic and transport prop- 
erties of these disordered lattices are investigated using ir-ir* 
orthogonal tight-binding (TB) model with nearest neighbor 
hopping 70 = 2.8 eV— and zero onsite energies. No varia- 
tion of the hopping elements with disorder is included in the 
model as bond-length variation does not exceed a few per- 
cent (cf. Tabled; all dependence on disorder stems from the 
ring statistics, which is expected to be the dominant effect 2 !. 
Fig. |TJd) shows the density of states (DOS) of the two disor- 
dered samples, together with the pristine case (dashed line) 
for comparison. Sample 1, which keeps 52% of hexagonal 
rings, displays several noticeable features, similar to those 
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FIG. 2: (color online) Normalized time-dependent diffusion coeffi- 
cients for two selected energies for both samples SI and S2. Inset: 
localization lengths as a function of the carrier energy 



which in the long time limit gives the conductivity through 
Eq. Q] Assuming an isotropic system in the x and y direc- 
tions, the 2D diffusion coefficient becomes D(t) = D x (t) + 
D y (t) = 2D x (t). All information about multiple scattering 
effects is contained in the time-dependence of D(t). 

Numerically, whatever the initial wavepacket features, 
D(t) starts increasing ballistically at short times, then 
reaches a maximum value which depends on the disorder 
strength, and finally decays as a result of quantum inter- 
ferences, the strength of which will dictate either a weak 
or a strong Anderson localization regime. The semiclas- 
sical quantities (elastic mean free path £ e (E) and semi- 
classical conductivity er sc ) are derived from the maximum 
of D{t) as i e {E) = D max (E)/2v(E) and a sc (E) = 
je 2 p(E)D max (E), respectively (with v(E) being the carrier 
velocity). 

The conductance of the amorphous stripes contacted 
to pristine graphene electrodes is computed using the 
Landauer-Biittiker approach^: 



found in previous studie s 15 ' 16 . First, the DOS at the charge 
neutrality point (Fermi level) is found to be increased by 
a large amount. Additionally, the electron-hole symme- 
try of the band structure is broken by the presence of odd- 
membered rings which generate quasibound states at reso- 
nant energies^ 6 .. The hole part of the spectrum is still remi- 
niscent of the graphene DOS, with a smoothened peak at the 
van Hove singularity, while for the electron part a second 
maximum appears close to the upper conduction band edge. 
By reducing further the ratio of even versus odd-membered 
rings (sample S2), the second maximum develops to a strong 
peak at about E = 2.5-jq while spectral weight at E = 370 
is suppressed. The redistribution of DOS at the upper con- 
duction band edge is a signature of odd-membered rings and 
its strength with increasing number of such rings relates the 
statistical distribution of rings with the DOS features. 

Transport Methodology.- To explore quantum transport 
in these topologically disordered graphene samples, a real- 
space order- A" quantum wavepacket evolution approach is 
employed to compute the Kubo-Greenwood conductivity^. 
Such method, pioneered in^t, has been successfully applied 
to many different types of systems, and in particular it has 
allowed to scrutinize Anderson localization in oxygen func- 
tionalized graphene—. The zero-frequency conductivity for 
carriers at energy E is computed as 



G(E) = G T(E) = ^-TV^t] 



(3) 



a dc = e 2 p(E) lira- : AX 2 (E,t) 



(1) 



where p(E) is the density of states and AX 2 (E,t) is the 
mean quadratic displacement of the wave packet at energy 
E and time t: 



AX 2 (E,t) 



Tr[5{E-U)\X(t)-X(0)f 
Tr[5(E - %)] 



(2) 



A key quantity in the analysis of the transport properties 
is the diffusion coefficient: D x (E,t) = f t AX 2 (E F ,t), 



where T(E) and t(E) are the transmission probability and 
transmission matrix, respectively, which can be computed 
from the Green's function G(E) in the contact region and 
the broadening T(E) of the states due to the interaction with 
the left and right electrodes. We calculate the conductance 
of the stripes, which are infinite and periodic in the direc- 
tion parallel to the interface with the graphene electrodes (y 
axis in Fig. 0I C ))- Despite the very large periodicity of our 
models, we perform a thorough sampling of the /c^-points in 
that direction^, to obtain the appropriate V-shaped conduc- 
tance of graphene in the thermodynamic limit. G is given 
per supercell of periodicity W=\ 1.4 nm. Note that, with this 
geometry, conductivity and conductance are related though 
a = ±G. 

Mean Free Path, Conductivity and Localization Effects.- 
Fig. |2] shows the time dependence of the normalized diffu- 
sion coefficient D(t)/ D max for two chosen energies and for 
both samples SI and S2. For energy E = —2jo, the short- 
time ballistic regime is followed by the saturation of the dif- 
fusion coefficient typically after 0.1 ps (for both samples SI 
and S2). From the saturation values, £ e (E) and a sc are de- 
duced and reported in Fig. [3] A striking feature is the very 
low value of the mean free path £ e , below 0.5 nm for the 
energy window around the Fermi level (where the DOS is 
considerably larger than that of pristine graphene). For neg- 
ative energies (holes) far from the charge neutrality point, a 
considerable increase of more than one order of magnitude 
in the mean free paths is observed. The increase occurs for 
smaller binding energies for sample S 1 than for sample S2, 
in good correlation with the changes observed in the DOS 
(which, around the van Hove singularity, deviates from the 
pristine graphene one more strongly for sample S2). 

We observe in the inset of Fig. [3] that a sc shows a min- 
imum value <7™ n about 4e 2 /irh, in agreement with the 
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FIG. 3: (color online) Elastic mean free path versus energy for the 
two samples. DOS of sample SI is also shown for comparison, in 
rescaled unts. Inset: semiclassical conductivity of corresponding 
lattices. 



theoretical limit in the diffusive regime, already confirmed 
for other types of disordef 3 -^. However, in contrast to 
prior studies, conductivity values near the minimum are ob- 
tained over an energy range of several eV around the charge 
neutrality point. This indicates that transport is strongly 
degraded in the amorphous network compared to pristine 
graphene, for which the conductivity increases rapidly when 
the Fermi level is shifted away from the Dirac point. The 
charge mobility, fi(E) = a sc (E)/en(E), with n(E) be- 
ing the carrier density, is found to be about 10 cm 2 V _1 s _1 
for n = 10 11 — 10 12 cm" 2 , orders of magnitudes lower 
than those usually measured in graphene samples 13 . Such 
low conductivity and mobility values should be measured at 
room temperature, where the semiclassical approximation is 
expected to hold. 

The obtained short mean free paths and minimum (semi- 
classical) conductivities indicate a further marked contribu- 
tion of quantum interferences turning the system to a weak 
and strong insulating system with temperature drop. Inter- 
ference effects are evidenced by the time-dependent decay 
of the diffusion coefficient D(t)/ D max as clearly seen in 
Fig. |2 although with large differences depending on the 
chosen energy. For E = —2^q, such decay is weak but 
clearly more pronounced for the S2 sample (which is more 
disordered than SI). In sharp contrast, localization effects 
are much stronger at the charge neutrality point and develop 
from much shorter timescale (few nanoseconds). These dif- 
ferences will be reflected in the corresponding localization 
lengths. 

Based on the scaling theory of localization, an estimate 
of the localization length of electronic states is inferred 
from£{E) = e e (E)exp(nha sc (E)/2e 2 'p: The results are 
shown in Fig. [2] (inset). The amorphous samples are con- 
firmed to be extremely poor conductors, with localization 
lengths as low as £ ~ 5 — 10 nm over a large energy window 



around the charge neutrality point. One also observes that £ 
can vary by more than one order of magnitude depending on 
the disordered topology of the sample and rings statistics. 




FIG. 4: (color online) Landauer-Biittiker conductance (for W= 1 1 .4 
nm) of two amorphous stripes contacted to graphene electrodes 
with L = 1.6 and 8.6 nm, respectively. The conductance of a 
pristine graphene contact with the same lateral size (11.4 nm) is 
shown for comparison. The inset shows the dependence of the con- 
ductivity on the stripe size L; symbols: calculated points; line: fit 
to <j{L) 
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To further confirm the localization lengths estimated us- 
ing scaling theory, we compute explicitly the conductance 
of the amorphous graphene stripes contacted with pristine 
graphene electrodes, using the geometry shown in Fig. [He), 
as a function of the length of the amorphous contact L. 
Fig. [4] shows the conductance curves for two stripes with 
L = 1.6 and 8.6 nm, respectively, compared to that of a 
graphene contact with the same lateral size in the supercell 
{W = 11.4 nm). It is clear that the conductance of the 
amorphous samples is greatly reduced with respect to that 
of graphene, and that the reduction is more pronounced as 
the length of the amorphous contact increases. Also, while 
the conductance for the stripe with the smallest length is 
relatively smooth, it becomes more noisy with increasing 
L. This reflects the transition from a diffusive system, in 
which the length of the amorphous contact is longer than the 
mean free path but shorter than the localization length, to 
a strongly localized one in which the localization length is 
shorter than the length of the amorphous region. 

From the variation of the Landauer-Biittiker conductance 
with L, we can extract reliable values of the localization 
lengths, as in the Anderson regime the conductance should 
decay as G(L) ~ e^ L ^. The inset in Fig. |4]shows the value 
of the conductivity (obtained from the Landauer-Biitikker 
conductance) for several stripes with L ranging from 1.6 to 
15.3 nm, averaged over an energy window of I.570 around 
the Fermi energy. A fit of the results to a(L) ~ w e ~ L ^ 
yields £ = 5.8 nm. This value is fully consistent with that 
obtained above using scaling theory for energies close to the 
Fermi level, and confirms that, in these amorphous struc- 
tures, strong localization effects should occur at low tem- 
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peratures at distances of less than 10 nm. Our results are 
consistent with the experimental ones from transport mea- 
surements by Zhou et al.—, which show values in the range 
between 0.1 and 10 nm for samples amorphizedby ion radi- 
ation. 

In conclusion, we have shown that topological disorder 
alone causes amorphous graphene to be a strong Anderson 
insulator. The increase of the density of states close to the 
charge neutrality point is associated with quantum interfer- 
ence which inhibits current flow at low temperatures. Very 



short mean free paths and localization lengths are predicted, 
in line with recent experimental evidence in graphene under 
heavy ion irradiation damaged. 
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